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ABSTRACT 

We present results from a 220-ks observation of the galactic supernova remnant (SNR) W49B 
using the Advanced CCD Imaging Spectrometer (ACIS) on board the Chanrda X-ray Observatory. 
We exploit these data to perform detailed spatially-resolved spectroscopic analyses across the SNR 
with the aim to investigate the thermodynamic properties and explosive origin of W49B. We find 
substantial variation in the electron temperature and absorbing column toward W49B, and we show 
that the mean metal abundances are consistent with the predicted yields in models of bipolar/jet- 
driven core-collapse SNe. Furthermore, we set strict upper limits on the X-ray luminosity of any 
undetected point sources, and we exclude the presence of a neutron star associated with W49B. We 
conclude that the morphological, spectral, and environmental characteristics of W49B are indicative 
of a bipolar Type Ib/Ic SN origin, making it the first of its kind to be discovered in the Milky Way. 
Subject headings: ISM: Supernova Remnants — X-rays: ISM — Supernovae: Individual (W49 B) — 
X-rays: ISM 



1. INTRODUCTION 



W49B (G43.3-0.2) is one of the most luminous Galac- 
tic super nova remnant (SNR) in X-rays (with Lx ~ 10 38 
Immler fc Kuntz||2005) and 7-rays, as observed 



erg s 



The physical origin of W49B's unusual traits remains 
uncertain. Two potent ial scenarios may account for the 
barrel shape of W49B flKeohane et al.||2QQ7| |Miceli et a~ 



by Fermi Large Area Telescope ( |Abdo et al.|2010|) . While 
radio maps show the SNR has a prominent shell ~ 4' in 
diameter, the X-ray emission is elongated in a centrally 
bright, iron rich bar which has two plumes at its east- 
ern and western edges (see Figure jl] |Hwang et al.j 2000 



Miceli et al]|2006| |Lopez et at 20(Md). The X-ray spec- 



trum ol W49B has strong emission lines from He- like and 
H-like ions of several metals with super solar abundances, 
indicative of their eject a-dominated origin. 

Previous studies of W49B have demonstrated the 
source has several unique morphological and spectral fea- 
tures which are not found in other Galactic SNRs. Fore- 
most, W49B is the only Galactic SNR wit h large-scale 
segreg a tion o f its nucleosynthetic products (Lopez et al. 



.20081 |Zhou et al.|2011| ): it could be due to a bipolar /jet- 
driven explosion ol a massive star, or it may be the 
result of an inhomogeneous ISM around a spherically- 
symmetric SN. The latter explanation is appealin g be- 
cause it applies to other sources (e.g., G292.0+1.8: |Park| 
et al.||2004| ). However, that scen ario would require un - 
realistic ejecta masses in W49B ( |Keohane et al. 2007). 
A bipolar /jet-driven Type Ib/Ic SIN origin is consis- 
tent with the median abundance ratios from a spatially- 
resolved sp ectroscopic analysis using XMM-Newton ob- 
servations (Lopez et al. 2009b), and it may account for 
the anisotropic distribution ol iron since heavy elements 
are preferentially ejected alon g the polar axis of the pro- 



2009b 2011): the iron synthesized in the explosion is . . . . , _ 6fY/ r . . _ 

miifg^ the Eastern half of the SNR, while the sudl events is estimated to be 0.7+ ; 6 % of local Type 



genitor in these explosions (|Khokhlov et al.||19 99; Maz- 
zali et al.|200T||Ramirez-Kuiz & MacFadyen|2010| |(J< 



et al.||201ip . Based on radio observations, the rate of 



intermediate-mass elements (like silicon and sulfur) are 
more spatially homogeneous. Secondly, the X-ray spec- 
trum is suggestive of rapid cooling uncommon in young 
SNRs, such as a strong radiative recombination contin- 
uum from iron (lOzawaet al.|2009| ). Finally, W49B is one 
of only two SNRs to date with detected X-ray emission 
l ines from the heavy elements chromium and manganese 



2010|), equating to ~1 per 



(Hwang et al. 2000 



Miceli et al. 



products ol explosive silicon burning . 



2006), nucleosynthetic 
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10 3 — 10 4 years per galaxy ( |Lopez et al.||2011[ ), with a 
subset of t hese explosions be ing hypernov ae or gamma - 
ray bursts fllzzard et al.|2004|[Podsiadlowski et al.|2004| ). 

In this paper, we investigate the explosive origin and 
thermodynamic properties of W49B based on a re- 
cent 220-ks Chandra X-ray Observatory Advanced CCD 
Imaging Spectrometer (ACIS) observation. In Section [2J 
we detail the observations and spatially-resolved spec- 
troscopic analyses of the data. In Section |3j we discuss 
the mult i wavelength evidence of the explosive origin of 
W49B and its implications. 

2. OBSERVATIONS AND SPECTROSCOPIC ANALYSES 

W49B was observed for 220 ks with Chandra ACIS 
on 18-22 August 2011 with the backside-illuminated S3 
chip in the Timed-Exposure Faint Mode (ObsIDs 13440 
and 13441). Reprocessed data were downloaded from 
the Chandra archive, and analyses were performed us- 
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Fig. 1. — Raw X-ray full-band (0.5-8.0 keV) exposure-corrected 
image of W49B as observed for 220 ks by Chandra ACIS. Point 
sources have not been removed to demonstrate the lack of any 
obvious pulsar or neutron star. Scale bar is V in length. 
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Fig. 2. — Integrated X-ray spectrum from the 220-ks Chandra 
observation taken in 2011 (red), with the integrated spectrum from 
previous Chandra data taken in 2000. The spectra have numerous 
prominent emission lines from He-like and H-like ions as well as 
a strong bremsstrahlung continuum. The high absorbing column 
toward W49B (~ 7 X lO^ 2 



cm 2 ) attenuates the soft X-rays below 



the Si xill line. In the 11-year interval between observations, the 
Si xill line flux decreased by ~18%, and the other emission lines 
have < 5% lower fluxes in 2011 compared to their 2000 values. 



ing the Chandra Interactive Analysis of Observations 
(ciao) Version 4.3. We extracted the integrated Chan- 
dra ACIS X-ray spectrum of W49B (Figure ^) using the 
CIAO command specextract for the central 2.9' around 
the SNR; a total of -10 6 net counts in the 0.5-8.0 keV 
band were recorded. The spectrum has a strong thermal 
bremsstrahlung continuum as well as prominent emission 
lines of silicon, sulfur, argon, calcium, chromium, man- 
ganese, and iron. 

We created exposure-corrected images of particular 
emission features by filtering the data to the narrow 



energy ranges corresponding to those features: Si xill 
and Si xiv (1.74-2.05 keV), S xv and S xvi (2.25-2.70 
keV), Ar xvn and Ar xvm (3.0-3.40 keV), Ca xix and 
Caxx (3.72-4.13 keV), and Fe xxv (6.4-6.9 keV). We 
continuum-subtracted these images by creating a contin- 
uum image (4.2-5.5 keV), scaling the image to reflect 
the continuum under the respective lines, smoothing all 
the images, and then subtracting the scaled continuum 
images from the line images. The resulting continuum- 
subtracted images are shown in Figure [3j Fe xxv is 
largely absent in the West of W49B, while the other ions 
are more homogeneously distributed. 

Chandra had previously observed W49B with ACIS-S3 
for 55 ks in 2000. In the 11-year interval between Chan- 
dra observations, the spectra have changed dramatically. 
In particular, the (unabsorbed) Si xill flux decreased 
by ~18%, and the other emission lines have < 5% lower 
fluxes in 2011 compared to their 2000 values (see Fig- 
ure^). Based on spatially-resolved spectral analysis, this 
line nux decrement is the greatest in the West of W49B, 
whereas the East has the smallest variation. Possible 
astrophysical explanations for this spectral evolution in- 
clude SNR expansion or changes in ionization state of 
the plasma; we defer further analysis and discussion of 
this spectral evolution to a future paper. 

2.1. Spectroscopic Analyses 

To examine the plasma temperature kT e and column 
density Ah across the SNR, we extracted spectra sys- 
tematically from__a grid of 713 regions of area 7.5" x 

extracted 
southeast of 

the SNR and were subtracted from the source spectra. 
Region sizes and locations were determined so that they 
have >400 net counts in the X-ray full band, the min- 
imum necessary to fit the bremsstrahlung shape accu- 
rately. Data were grouped to a minimum of 10 counts 
per bin, and sp ectra were fit using XSPEC Version 12.7.0 
( Arnaud|1996 ). We modeled the spectra as an absorbed, 
variable abundance plasma in collisional ionization equi- 
librium (CIE) using the XSPEC c omponents phabs and 



7.5" (see Figure [4]). Background spectra were 
from a source- free, 2' diameter region so 



vmekal (Mewe et al. 1985, 1986 Liedahl et al. 1995) 



The absorbing column density Ah, electron temperature 
kT e , abundances (of Si, S, Ar, Ca, and Fe, with solar 



abundances from Asplund et al.|2009| ), and normalization 
were allowed to vary freely during the fitting. Although 
previous studies employed t wo CIE plasmas in model 



Miceli et al. 2006 Lopez 



le use ol two components 



ing the X-ra y spectra (e.g., 
et al.||2009b| ), we find that t 
instead of one for our small regions does not statistically 
significantly improve the reduced chi-squared values of 
our fits. Therefore, we employ a one-temperature CIE 
plasma model to fit the spectra of our regions. 

The results of the spatially-resolved spectral analyses 
of the small, 713 regions are presented in Figure [4j where 
we map the best-fit values of Nn and of kT e . The fits 
were statistically good: the mean reduced chi-squared 
X 2 value was 1.11 with 87 degrees of freedom, with a 
1-cr range of \ 2 =1.02-1.19 with 53-101 degrees of free- 
dom. The temperature and absorbing column density 
vary substantially across the SNR, with values ranging 
N u -4-12 xlO 22 cm" 2 and kT e ~ 0.7-2.5 keV (with esti- 
mated uncertainties of 10% in both parameters for each 
region). The A^h is enhanced in the northeast and south- 
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Fig. 3. — Exposure-corrected, continuum-subtracted images of Si xm and Si xiv, S xv and S xvi, Ar xvn and Ar xvm, Ca xix and Ca 
XX, and Fe xxv. Images have been smoothed with a Gaussian kernel of a =5 pixels. The scale bar is V in length; North is up and East is 
left. Fe xxv emission is weak in the West of W49B, while the other ions are more homogeneously distributed. The "hole" in the Si image 
in the Northeast of the SNR is caused by the large absorbing column JVh ~10 23 cm -2 there. 
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Fi g. 4. — Left: Grid of 713, 7.5" x 7.5" regions where spectra were extracted. Contours show the location of warm H2 (at 2.12 /im; 
from Keohane et al. 2007J. In the East, the ejecta are impacting the molecular material, while the ejecta are likely freely expanding in 
the West. Maps ot best-lit values from spectral analyses for absorbing column density N-r (middle) and electron temperature kT e (right). 
Overplotted are the X-ray full-band image contours to guide the eye. 



west of the SNR, with many regions having Nn > 8x 10 22 
™-.-2 This result is consistent with the work of Lacey 



et al. (2001) who showed there may be intervening ab 
sorbing material along the line of sight toward W49B 
in the southwest using low- frequency radio observations. 
Overall, our measured Nn valu es are greater than those 
found in previous studies (e.g., |Hwang et al" 2000 found 
7V H ~ 5.0 x 10 22 cm -2 using integrated ASCA spectra) 
be cause we have emplo yed the revised solar abundances 
of |Asplund et al.||2009| We find a temperature gradient 
across W49B, with elevated kT e toward the Eastern side 
of the SNR. The plasma is likely hotter there because it 
is impacting nearby molecula r material to the Eas t (see 
the contours in Figure |4j left; Keohane et al.||2007[ ). 

To assess the abundances of the ejecta metals, we repli- 
cated the spatially-resolved spectroscopic analyses de- 
scribed above (i.e., modeled the data as an absorbed, 
single CIE plasma with variable abundances) using 136 
regions of area 15" x 15". Regions were selected such that 
they had > 2500 net full-band counts, the minimum nec- 
essary to accurately model the line emission, and the 
mean number of full-band counts per region was 5500. 
We set the lower limit on the count rate by produc- 
ing confidence contours of the abundances for regions 
of varying count rates, and we found that regions with 
<2500 net full-band counts did not yield reliable con- 
straints. We opted to use 15"xl5" regions larger area 
regions were not adequately fit with a single-temperature 
plasma (i.e., larger regions may superpose distinct plas- 
mas of different temperatures). The mean reduced chi- 
squared of the fits was x 2 = 1-50 with 161 degrees of 
freedom, with a 1-a dispersion range of x 2 =1.29-1.71 



with 92-232 degrees of freedom. 

In Figure [5j we show histograms of the abundances of 
Si, S, Ar, Ca, and Fe given by our best fits of the 136 
regions. The mean abundances of these elements rela- 
tive to solar by number and their 1-a range are 3.0^"J, 



2.7±1.1, 2.9 



+0.9 



4.1 



+1.5 



+4.8 



14 , and 9.1_ 5 



respectively. We 

estimate the uncertainties of the abundances in the fits of 
individual regions to be ~20-30%. Therefore, we expect 
that the large dispersion in abundances reflects the non- 
uniform distribution of the elements across the SNR. In 
particular, the large variance in the Fe abundance arises 
from the fact that the iron is absent in the Western half 
of W49B, where many regions have Fe at or below their 
solar values. As the temperatures (from the map in Fig- 
ure [4| in the West are high enough to produce prominent 
emission from Fe x xv if iron is present (see Figure 23 of 
Lopez et al.||2009b ), this result indicates that iron is de- 
pleted in the West. 

3. DISCUSSION OF THE EXPLOSIVE ORIGIN OF W49B 

In Figure [6j we plot the mean abundances by mass 
and their 1-a range relative to iron of silicon (Si/Fe ~ 
), sulfur (S/Fe « 0.10±g;g|), argon (Ar/Fe « 



0.22™), sulfur 

7+0.012 



0.027lo.oi6)> and calcium (Ca/Fe « 0.0341^^), from 
the fits of the 136 regions. For comparison, we show the 
abundance ratios predicted by six core-c ollapse models: 
four e nergetic, spherical explosions (from [Nomoto et aTj 
2006) an d two aspherical bipolar explosions (25 A and 
25B from |Maeda fc Nomoto| |2003). Generally, asymme- 
tries in the explosion alter the chemical yield of bipolar 
explosions relative to typical CC SNe: while typical CC 
events eject ~ 0.07 — O.15M of 56 Ni, bipolar explosions 
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Fig. 5. — Histograms of the abundance ratios by mass obtained from the spatially-resolved spectroscopic analyses across W49B. 
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have much greater nickel yields which increase with as- 



phericity, progen itor mass, and explosion energy (Umeda 
& Nomoto 2008). The primary difference between the 



two aspherical models in Figure [6| is their jet opening an- 
gles, 0: 25 A has = 15° and 25Bhas = 45°. For a fixed 
kinetic energy, larger means less mechanical energy is 
released per unit solid angle, affecting the luminosity of 
the explosion and the nucleosynthetic products. 

The ratios we obtain for W49B agree better with the 
predictions of the aspherical models than those of the 
spherical models. The Si/Fe of the aspherical model with 
larger opening angle, = 45°, is within ~5<r range of 
our fit values, while the spherical explosions have Si/Fe 

> 17cr away from the measured ratios. The S/Fe fit values 
are also within 2a of the prediction for the aspherical 
model with = 45° and within 4.5a range of that for the 
= 15° aspherical model, whereas the measured S/Fe is 

> 10a from the ratios of the spherical explosions. For 
Ar/Fe, the 6 = 45° model gives a value roughly equal 
to our mean estimated Ar/Fe, and the 6 = 15° model 
is ~1.5<j from the obtained value. By comparison, the 
spherical models are > 5<r from the measured Ar/Fe in 
W49B. Finally, both aspherical models give Ca/Fe ratios 
within the la range of those estimated in W49B, and the 
spherical models give values > la above the measured 
means. Based on the results for each of these four ratios, 
we find that a bipolar origin of W49B is most consistent 
with the abundances found across the SNR. 

To exemplify an SNR that has abundances consistent 
with spherical model predictions, we measured Si/Fe 
and S/Fe in several regions of G292.0+1.8, an SNR 
thought to be from a sphe rically-symmetr ic CC explosion 
within a complex CSM ( |Park et aT]|2004[ ). We extracted 
spectra from 26 regions id entified as ejecta substruc- 



tures in [LopezetaL (2011) and modeled these data as 
an absorbed, variable abundance plane-parallel shocked 
plasma with constant temperature (with XSPEC model 
components phabs and vpshock). We found mean abun- 
dance ratios of Si/Fe « 1.54 and S/Fe w 0.62. These val- 
ues are consistent with the M = 25 M ft a nd M = 30 M 
spherical models of |Nomoto et al.|( 2006). 

The recent 220-ks Chandra AClfe observation has fa- 
cilitated a detailed study of the ejecta and thermody- 
namic properties across W49B, and the measured abun- 
dance ratios are close to those predicted for a bipolar /jet- 
driven explosion of a 25 M progenitor. By contrast, 
these ratios are inconsistent with the values for symmet- 
ric CC SNe of different explosion energies and progenitor 
masses. Although the bipolar SN models do not match 
the measured abundance ratios precisely, the available 
models of aspherical CC SNe are limited, and the work 
presented above demonstrates the need for more simula- 
tions exploring the effects of asphericity, progenitor mass, 
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Fig. 6. — Plot of mean abundances (by mass) and their la range 
of values relative to iron of silicon, sulfur, argon, and calcium. For 
comparison, we show the abunda nce ratios predicted b y models 
of four spherical explosion s (from Nomoto et al. 2006 ) and two 
aspherical explosions (from Maeda & INomoto 2UU3 ). 

and explosion energy on nucleosynthetic yields. 

In addition to its anomalous abundances, W49B has 
other observational signatures consistent with a bipolar 
origin. The X-ray line emission of W49B is the most 



elliptical and elongat ed when com 
X-ray bright SNRs (|Lopez et al. 



pared t o other young, 
2009a), and the bar- 



rel shape of W49B is what would be expected from a 
bipolar SNR if the jet was oriented near the plane of 
the sky. Additionally, the Fe xxv is morphologically dis- 
tinct from the lower- Z ions: the iron is more asymmetric 
and is segregated from the lighter elements ( Lopez et al. 
2009b|). Furthermore, as the temperatures in the loca- 



tions without prominent Fe xxv emission are sufficient 
to produce that feature, the anomalous distributio n of Fe 
xxv likely a rises from anisotropic ejection of iron (Lopez 
et al.|2009b| ) . The morphological differences of iron from 
lower- Z elements is consistent with a bipolar origin: the 
enhanced kinetic energy at the polar axis causes nickel to 
be more efficiently synthesized there, while perpendicular 
to the rotation axis, mechanical energy decreases, caus- 
ing lo wer- Z elements to al so be ejected more isotropically 
(e.g., [Mazzali et al.|[2QQ5] ) 



The environment ol W49B is also suggestive of a mas- 
sive progenitor sta r that would be n ecessary for a Type Ic 
bipolar explosion. |Keohane et al.| ( |2007[ ) imaged W49B 
at near-infrared wavelengths and found a barrel-shaped 
structure in the 1.64 jira [Fe n] which they interpret 
as coaxial rings from wind-blown material of a massive 
star. These authors also found warm, shocked H2 at 
2.12/im enclosing the SNR to the East and Southwest, 
while cooler 13 CO gas has been detected at the North 
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Fig. 7. — Plot of total photon luminosity at infinity L 

oo versus 

age for nearby neutron stars and predictions from models of neu- 
tron star cooling. The light blue diamond is our upper limit for a 
neutron star in W49B; dark blue points are detect ions and red tri- 
angles are upper limits, from l T suruta et al. (2009). Shaded curves 
are model predictions of Page et al. ( 2UU6fior neutron star cooling 
with a heavy element envelope (dark gray), a maximum amount of 
light elements (medium gray), and an intermediate amount of light 
elements (light gray). Our upper limit for W49B (logLoo < 31.6) is 
two to three orders of magnitude below the predictions for a 1000 
year old source. 



and West of W49B ( [Simon et aLl|2QQl[ ). This molecular 
gas is likely from the cloud where the massive progenitor 
of W49B formed. Moreover, W49B i s one of the mos t 
luminous SNR in GeV gamma rays dAbdo et al.pOlOL 
and the brightest SNRs at these energies are interacting 
with molecular clouds because pion decay of relativistic 
ions is most efficient when they are accelerated in shocks 
advancing through dense material. 

Finally, as the remnant of a bipolar/jet-driven explo- 
sion, W49B might be expected to have formed a black 
hole or a magnetar, and thus no neutron star or pul- 
sar should be detected in the source. From our 220-ks 
Chandra observation, we estimate any undetected neu- 
tron stars would have an upper limit X-ray luminosity 
of L x < 3.7 x 10 31 erg s" 1 in the 0.1-10 keV band (as- 
suming a 0.5 keV blackbody). Figure [7] compares our 
upper limit in W49B to current detections and u pper 
limits of therma l emission in nearby neutron stars (Ka- 
plan et al. 2004j Tsuruta et al. 2009). Furthermor e, we 
plot model predictions ol neutron star cooling (from Page 



et al. 



2006). For a 1000 year old point source, the Page 



models predict total photon luminosities at infinity of 
logLoo ~ 33.3 — 35. Thus, if W49B hosts a neutron star, 
our limit requires the central source to be less luminous 
than those currently known and 2-4 orders of magnitude 
lower than pre dictions of neut r on sta r cooling models. 

Furthermore, Gorham et al. ( 1996 ) searched for pulsa- 
tions in W49B using Arecibo, and they found no pulsa- 
tions with periods >0.5 ms down to flux densities of 1.2 
mJy at 430 MHz across the entire SNR. Finally, the spec- 
trum of the bright gamma-ray emission detected with 
Fermi-LAT fr om W49B is inc onsistent with it arising 
from a pulsar ( Abdo et al.||2010[ ). Given the strict upper 
limits on the presence ol a neutron star /pulsar, it is rea- 



sonable to assume that the compact object produced in 
the explosion might be a black hole. 

Therefore, the morphological, spectral and environ- 
mental evidence supports a bipolar/jet-driven SN origin 
of W49B, and it is the first SNR of this kind to be identi- 
fied in the Milky Way. Given the estimated rate o f these 
explosions is ~1 per 10 3 — 10 4 years per galaxy (Lopez 
et al.pOlT ), it is feasible th at a few of th e 274 currently 
known SNRs in the Galaxy ( Green|[2009 ) originate from 
this kind of SN. 

In the future, this result can be verified two ways. 
If W49B was a bipolar explosion, it should have copi- 
ous intermediate-mass elements produced through hy- 
drostatic burning of a massive star (like O and Mg). 
While the high Nn toward W49B precludes detection 
of emission lines from these metals in the optical and X- 
ray, we predict emission features from elements should 
be detectable at near - infrared wavelengths ( as observed 
in e.g., Cassiopeia A: Gerardy fc Fesen||2001 ). Secondly, 
if W49B was jet driven, the metals in the central bar 
should be moving at velocities greater than the more 
homogeneously distributed, lower- Z plasma. Therefore, 
using spatially-res olved, high-resolution X-ray spectrom- 
eters like Astro-H ( |Takahashi et al.|2010 ), one could mea- 
sure the Doppler shifts ol the jet plasma and verify its 
higher speeds than other eject a material. 
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